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Abstract 

In this article we discuss the well-posedness of a mathematical model 
that is used in the literature for the simulation of Lithium-ion (Li-ion) 
batteries. First, a mathematical model based on a macro-homogeneous 
approach is presented, following previous works. Then it is showed, from 
a physical and a mathematical point of view, that a boundary condition 
widely used in the literature is not correct. Although these errors could 
be just sign typos (that can be explained as carelessness over d/da; versus 
d/dn, with n the outward unit vector) and authors using this model prob¬ 
ably use the correct boundary condition when they solve it in order to do 
simulations, readers should be aware of the right choice. Therefore, the 
deduction of the correct boundary condition and a mathematical study of 
the well-posedness of the corresponding problem is carried out here. 

Keywords: Lithium-ion batteries; mathematical model; partial differential 
equations, boundary conditions; uniqueness of solution; existence of solution 


1 Introduction 

Lithium-ion (Li-ion) batteries have become very popular in the last years as a 
source of energy in multiple portable electronic devices. A mathematical model 
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showing the key factors of the battery operation can be very helpful for the 
design and optimization of new models and also for the real time control of its 
performance. 

Based on a macro-homogeneous approach developed by Newman (see 0) 
several mathematical models have been developed for these purposes (see 01113 
0[TTl|Tl[Tni[Tlini[Z]) which include the main physics present in charge/discharge 
processes. A fully mathematical model is presented in Sec.0 including a system 
of boundary value problems for the conservation of Lithium and conservation of 
charge in the solid and electrolyte phases, together with an initial value problem 
for the conservation of energy. In the literature (see [T], [B], 0, [TU], [H], [H]) 
one can find numerical computations of the model (or simplifications of it), with 
parameters corresponding to actual devices, that help to highlight the structure 
of this highly coupled model and show the relevance to the applications. 

Over the last years, several authors have written articles in different journals 
including a boundary condition that is not correct (see Remark^) . Although the 
authors probably use the correct boundary condition when solving the model 
in order to do simulations, the reader should be aware of the right choice. In 
this article it is shown why that condition is not only physically incorrect (see 
Sec. [3]) but also mathematically incorrect, since it is proved (see Sec. 13 Remark 
El) that the corresponding system of boundary value problems modelling the 
conservation of charge does not have any solution (therefore the system is not 
well-posed). Some results regarding the uniqueness and existence of solution 
(with the correct boundary condition) for a simplified version of the model are 
shown in Sec. |3 

2 Mathematical model 

2.1 Generalities 

A typical Li-ion battery cell has three regions: A porous negative electrode, 
a porous positive electrode and an electron-blocking separator. Furthermore, 
the cell contains an electrolyte, which is a concentrated solution containing 
charge species that can move all along the cell in response to an electrochemical 
potential gradient. 

The negative electrode is an intercalated lithium compound usually made 
from Carbon (typically graphite), with Li^Cg active material. Here y € [0,1] 
is the stoichiometry value of the material, which changes during charge and 
discharge. For instance, during discharge, lithium ions inside of solid Li^Cg 
particles diffuse to the surface where they react and transfer from the solid 
phase into the electrolyte phase (see [T3 El])- During charge they follow the 
opposite way. Theoretically, in a fully charged Li-ion battery this compound is 
LiCg (lithiated graphite; Li saturation; y = 1), in a semi charged/discharged 
battery it is LiyCg, with y G (0,1) and in a fully discharged battery it is just 
carbon (Li depletion; y = 0). In practical operating cases y never attains the 
extreme values y = 0 or y = 1. The corresponding negative reaction equation 
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is the following: 


yLi++ye + 6C <—S> Li^Ce (0 < y < 1). 


The positive electrode is usually a metal oxide or a blend of multiple metal 
oxides (see m) such as lithium cobalt oxide (Lii_yCo 02 ), lithium iron phosphate 
(Lii_yFeP 04 ), or lithium manganese oxide (Lii_yMn 2 04 ), with y € [0,1]. For 
instance, during discharge the positively charged ions travel via diffusion and 
migration through the electrolyte solution to the positive electrode where they 
react and insert into solid metal oxide particles (see, e.g., [nmu). During 
charge they follow the opposite way. Theoretically, in a fully charged Li-ion 
battery y = 1 (Li depletion), in a semi charged/discharged battery y G (0,1) and 
in a fully discharged battery y = 0 (Li saturation). Again, in practical operating 
cases y never attains the extreme values y = 0 or y = 1. The corresponding 
positive reaction equation for the examples showed above are the following: 


LiCo02 ^—5- Lii_yCo02 -I- yLi"^ -|- ye (0 < y < 1), 

LiFeP04 t—S’ Lii_yFeP04 -f yLi^ -I- ye“ (0 < y < 1), 
LiMn204 ■<—> Lii_yMn204 -I- yLi’*’ -|- ye~ (0 < y < 1). 


Therefore, considering both, negative and positive reaction, the corresponding 
total reaction equations are the following: 


LiCo02 -f 6C t^ Lii_j^Co02 -f Lij^Ce (0 < y < 1), 

LiFeP 04 + 6C t^ Lii_yFeP 04 -f Li^Ce (0 < y < 1), 
LiMn 204 -I- 6C t — )■ Lii_yMn 204 -I- Li^Ce (0 < y < 1). 


A ID electrochemical model is considered for the evolution of the Li con¬ 
centration Ce{x,t) (mol / m^) and the electric potential (V) in the elec¬ 

trolyte and the electric potential cj)s{x,t) (V) in the solid electrodes, along the 
^-direction, with x G (0, L) and L = Li + S + L 2 being the cell width (m). We 
assume that (0,Li) corresponds to the negative electrode, (Li,Li -I- S) corre¬ 
sponds to the separator and (Li -|- 5, L) corresponds to the positive electrode. 
This is coupled with a ID microscopic solid diffusion model for the evolution 
of the Li concentration Cs{x;r,t) in a generic solid spherical electrode particle 
(situated at point x G (0, Li) U (Li -I- 6,L)) along the radial r-direction, with 
r G [0,i?s] and Rs (m) the average radius of a generic solid active material 
particle. This ID approximation is valid since the characteristic length scale 
of a typical Li-ion cell along the x-axis is on the order of 100 ym, whereas the 
characteristic length scale for the remaining two axes is on the order of 100,000 
ym or more (see DP)- Rs can be different at each electrode and therefore we 
consider 
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The percentage of available local energy at time t (s) and radius r of a generic 
negative electrode particle situated at point x G (0,Li) in the cell x-direction 
is the the same as its stoichiometry value y = y{x, r, t), which can be computed 
as 


y{x,r,t) 


Csjx] r,t) 

^s, —,max 


where Cs,-,max (mol m“^) is the maximum possible concentration in the solid 
negative electrode and Cs{x; r, t) (mol m“^) is the the Li concentration at time t, 
radius r and point x. Therefore, the bulk state of charge (SOC) for the negative 
electrode (it can be also done for the positive electrode but both are related 
and, therefore, it suffices to use only one of them) is 


SOC(t) 







( 1 ) 


We point out that SOC is an nondimensional quantity that can be used as an 
indicator of the available energy in the cell. Measuring Cs{x] r, t) is not easy, but 
we can use the mathematical model below (see, e.g., [U [2l [TTl [12] ) to compute 
it. From a theoretical point of view SOC could be a value between 0 and 1 but, 
as mentioned above, in practical operating cases it never attains the extreme 
values 0 or 1. Therefore intermediate values ?/o%i -2^0%; 2/ioo%! ■ 2 ^ 100 % ^-re used for 
each electrode to refer to 0 % or 100 % SOC. 


2.2 A full mathematical model 

Based on the models appearing in the literature (see [3ll[5lll[TI||l2l[Tn|[Tl 
lEl) and assuming constant diffusion and activity electrolyte coefficients, a 
full mathematical model for the performance of a battery, also including heat 
transfer dynamics, is given by system of equations (ED-®: 


9ce d 
dt ^dx 


, dCe 

' dx 


1 — ^0 

—in (0,L) X (0,tend), 


f)c ric 

^(0,t) = ^(L,t) = 0, te(o,w), 

Ce(a;, 0) = Ce,o(a;), a;e(0,L), 


( 2 ) 


For each x G (0, Li) U (Li + 5,L) : 




d f ^dc, 
r 


dcs 
~dt 

dcs 

—T ^ a 

or or 

_ Csix;r,0) = Cs,oix;r), 


^ = 0, in (0,i?s) X (0,tend), 


dr dr 
{x;0,t)=0, -Ds^{x;Rs,t) = 


dcs, „ Rb{x) .Li 


(3) 


iSs{x)F 


j\ tG(0,tend), 
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For each t G (0, tend) '■ 


-i ("•'‘s'"w) w 

^(0,,) = ^(L,,) = 0, 

For each t G (0,<end) : 

in (0, Fi) U (Fi+5, L), 

e.(0)a(0)^(0,t) = e^{L)a(L)^{L,t) = 

^(ii,t) = ^(ii + ^,i)=0, 


j ~ — Tamb^ + 9r + 9j + 9c + 9e, t G (0, tend), 

1 r(0) = To, 

In the above system of equations Ce = Ce{x,t) (mol / m^) at time t (s), 
Cs = Csix;r,t) (mol / m^), (/)e = (t>e{x,t) (V), (j)^ = (t)s{x,t) (V), T = r(t) is the 
temperature of the cell (K), I = I (t) is the applied current (A), 

{ £e,- ifa;G(0,Li), 

^e,sep if X G (Fi, Li S), 

£e,+ if a; G (Ai + S, L) 

is the volume fraction of the electrolyte, p is the Bruggeman porosity exponent 
(nondimensional constant), De is the electrolyte diffusion coefficient (m^ s“^), 
(dimensionless and assumed here to be constant) is the transference number 
of Li+, 


De = De{x) = 


De- ifxG(0,Ai), 


" ' \ De,+ ifxG(Li+5,L) 

is the solid phase Li diffusion coefficient (m^ s“^), 


£s = £s(a;) = 


£s- ifa;G(0,Li), 
£sj+ if x G (Li + 5, L) 


is the volume fraction of the active materials in the electrodes, k = k (ce(x, t),T(t)) 
is the electrolyte phase ionic conductivity (S m“^), A (m^) is the cross-sectional 
area (also the electrode plate area), 


CTs = aeix) = 


as- ifxG(0, Li), 

crs,+ if a; G (Li -I- 6, L) 
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is the electrical conductivity of solid active materials in an electrode (S m 
and 

jhi ^ (l)six, t), (l)e{x, t),Cs{x; Rsix), t), Ce(x, t), 

is the reaction current resulting in production or consumption of Li (A m“^) 
at point X and time t. For the Butler-Volmer equation is usually used (see, 
e.g, [llSlIIIllIlIinilllll]): 





3£s(a:) . 

i?s(x) 


exp 


a.F 


-V 


exp 


-OLc.F 


-rj 




RT 'J " \ RT 
iixe (0,Li)U(Li-Hd,L), 


(7) 


, 0 if X G (Li, Li -I- S) 


(here, for the sake of simplicity, we have considered the solid/electrolyte inter¬ 
facial film resistance to be zero and therefore is not included in the above equa¬ 
tion), where as(x) = is the specific interfacial area of electrodes (m“^), 

^0 — ^0 (x^ Cg, Ce), 


Zo(x, Cg, Ce) 


fc_(Ce)“'‘(Cg _,niax - Cg)“'‘(Cg)“'= ifxG(0,Li), 

fc+(Ce)“'*(Cs, + ,max - Cs)““(Cs)“= if X G (Li + 6, L) 


is the exchange current density of an electrode reaction (A m“^), fc_,fc-|_ are 
kinetic rate constants (A mol“^“'*““‘=), Oa, otc (dimensionless con¬ 

stants) are anodic and cathodic transfer coefficients for an electrode reaction, 

■q = q(^x,4>six,t),(j)eix,t),Cs{x;Rs{x),t),T(t)^, 


•n 




(^g - (()e - U{x,Cs,T), if X G (0,Li) U (Li + S,L), 
0 if X G {Li,Li + 5) 


is the surface overpotential (V) of an electrode reaction and U{x,Cs,T) is the 
equilibrium potential (V) at the solid/electrolyte interface (i.e. the open circuit 
voltage - OCV). A way of expressing U is (see [5]): 


U{x,c,T) = 


U_ 

U+ 


Cs, —,max 
C 


+ 


dU_ 

~dT 

dU+ 


Cs, + ,max / dT Vcs.+. 


-)(T-T,ef) ifxG(0,Li), 
-)(T-T,,f) ifxG (Li + d,L), 

c/ ^ ^ 


dU- dU- 

where U-,U+, are functions tipically obtained from fitting experi¬ 

mental data and Cg^+^max (mol m“^) is the maximum possible concentration in 
the solid positive electrode. 
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Remark 1 For low overpotential cases Butler- Volmer equation can be sim¬ 
plified to the following linearized version (see m): 




3 5 06 5 Cs5 ^6 5 - 


T = 


SSsix) . (aa + ac)^' ^/n r ^M,'r , x r^ 

p . X »o- ^ - V, tfxG{0,Li)U{Li+S,L), 

Rs{x) R 1 

0, if X G (Li, Li + (5). 


Regarding heat transfer Eq. dh]), M (kg) is the mass of the battery, Cp (J 
kg“^ K“^) is the specific heat capacity, h (W K“^) is the heat transfer 
coefficient for convection, Ag (m^) is the cell surface area exposed to the convec¬ 
tive cooling medium (typically air), Tamb is the (ambient) temperature of the 
cooling medium, 

= qr{t)= A / j^^pdx 

Jo 

is the total reaction heat, gj = qj{t), 


gj(t) = Ass-a- 


dx 


2 i-L 

dx -|- Aes,+a+ 






R T 


dfs 

dx 

91nce 

dx 


dx 


dfe 

dx 


dx 


is the ohmic heat due to the current carried in each phase and the limited 
conductivity of that phase. 


Qc = qdt) = I{tf^ 

is the ohmic heat generated in the cell due to contact resistance between current 
collectors and electrodes, R{ {ft m^) is the film resistance of the electrodes and 


Qe = qe{t) 



■.j dU- 

dT 

dU+ I 
dT \ 


Cs(x, Rg , t) 


Cs. — 


Cg (X, Rg | , t) 
Cs,-t-,niax 


j dx 
dx 


is the reversible heat caused by the reaction entropy change (see m and [5] for 
a simpler formulation of qe(t)). 

Heat sources gr, gj and gc are always positive and, as explained in [11], the 
second term inside the last integral of gj is generally negative. On the other 
hand, ge can be either positive or negative. 

Remark 2 The second term on the left hand side of system w is often written 
in the literature using (1 — t+) (see Refs, and or 2(1 —1+) (see Refs. [3 
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[m [n [1210 instead of {1 — which is the correct term (see pag. 250 
of Ref. We remind that we are assuming here that t^ is constant and it is 

given (see m) by 

0 ^ Do,+ 

+ L»o,-+^o.+ ’ 

with Z)o,-i^o,+ transport coefficients related to the diffusion of Lt^ cations and 
the corresponding anions, respectively. 

Remark 3 Many authors (see, e.g., m and m) use some physiochemical pa¬ 
rameters involved in the model depending on T and use the Arrhenius equation 
to express such a dependency, which needs the corresponding activation ener¬ 
gies. The dependency of some parameters on T can sometimes be estimated. For 
instance in Ref. m the authors give formulae for K{ce,T), De{c,T) andt^{c,T). 

Remark 4 This system of equations does not have uniqueness of solution (if 
4>s{x,t) and 4>e{x,t) are solutions then <f>s{x,f) + c{t) and 4>e{x,f) + c{f) are 
also solutions, for any function c(t)). A way of avoiding that is to set a refer¬ 
ence value of 4>s{x,f) or (l)e{x,t) at some point x. For instance we can impose 
fsiOjt) — 0 for any t € [0,tend] (see Remark\T^. Some results regarding the 
existence and uniqueness of solution are showed in Sec. 0 

After solving the above model, we can get SOC(t) by using (0) and the cell 
voltage V{t) (A), at time t, given by 

V{t) = UL,t)-UO,t)-^nt)- ( 9 ) 

The battery pack voltage is calculated by multiplying the single cell voltage of 
Eq. (|ni) by the number of serially connected cells in the battery. 


3 Correct boundary conditions in ([5]) 


Remark 5 In ]11[ \12\. \1(K 0/, authors use the following boundary conditions at 
X = 0 and x = L, instead of those presented in m: 


d(t>s, , dcfs, , lit) 


( 10 ) 


In 0/ authors use the following boundary conditions at x = 0 and x = L: 




(11) 


In Sec. ro it is proved that both conditions are not correct from a physical 
point of view, because the resulting system of equations does not model a battery 
performance. Then, in Sec. \f.l\ it is proved that they are also mathematically 
incorrect, because, in fact, the corresponding system of boundary value problems 
modelling the conservation of charge does not have any solution (therefore the 
system is not well-posed) unless I{t) = 0 (see Remark^. 
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3.1 Deduction of system (1^ 

Given an external current lit) (A) applied to the battery (/(t) > 0 when the 
battery is discharging), according to Kirchoff’s law I = Is + Ig, where Ig, Ig are 
the current (A) in the solid electrode and in the electrolyte, respectively. Now, 
by Ohm’s law, 

Is{x,t) = -Aa{x)^^{x,t) = I{t) - Ig{x,t)- (12) 

Then, for each t G (0,tend), taking into account the porosity nature of the 
solid electrodes (only a fraction of its volume contribute to its conductivity) and 
using that 

r)T 

(13) 

the following equation for the conservation of charge in the electrode solid phase 
is obtained: 


= in (0,Li) U (Li + (5,L). 

We point out that Is{x, t) = 0 (and therefore Ig{x, t) = I[t)) for all x S (Li, Li + 
S) and lg{0,t) = Ie{L,t) = 0- This provides, using (fT^ . the following boundary 
conditions: 

s.(0)a(0)^(0,t) = eg{L)aiL)^{L,t) = 

-^^(Ti,t) — "^^(Ti + o,t) — 0; 

which completes the deduction of system ©■ Therefore, if for instance the 
battery is discharging, > 0 (respectively, < 0) if x € (0, Li) (respectively, 
if X G (Li + S,L)) and the graph of foi' some time t G (t,tend) could be 

qualitatively (not necessarily quantitatively) similar to that showed in Figure [TJ 


3.2 Deduction of system (]4l) 

Neglecting deviations of the electrolyte solution from ideal behavior 

Ie{x,t) = -AK^^{x,t) A{1 -2t^)^K-^\n (ce(x,i)) 

(see 0). The first term on the right hand side is due to Ohm’s law (as in (fT^ l 
and the second term accounts for concentration variations. 

Then, for each t G (0,<end), taking into account the porosity nature of the 
solid electrodes and using again 0, the following equation for the conservation 
of charge in the electrode solid phase is obtained: 


A 

dx 



+ ( 1 - 24 ) 


RT d 



= J 


Li 
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in (0,L), 
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_1 -^^^^- 

0 0.2 0.4 0.6 0.8 1 

x/L 

Figure 1: Typical Graph of for some time t G (t,tend). 


together with boundary conditions /e(0,t) = Ie{L,t) = 0. i.e., 

-K(Ce(0,t))^(0,t) + (1 - 2<° )^K(Ce(0,t))^ln(ce(0,t)) = 0, 

-K{CeiL,t))^{L,t) + (1 - 2t°^_)^K{CeiL,t))-^ln {Ce{L,t)) = 0, 
which, taking into account the boundary condtions in ([2]), are equivalent to 




4 About existence and uniqueness of solution of 
system (jl])—([5]) 

Let us study some properties regarding the existence and uniqueness of solution 
of (|4l)-([5|). We will only focus here on the equations involving conservation of 
charge. 

Even if coupling with the concentrations of Lithium is not included in the 
analytic work, the system under studied is interesting and not trivial at all. It 
is interesting because suitable numerical schemes of the full model will probably 
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need to solve this kind of systems at each time step of a time discretized version 
of the full model. It is non trivial because it is a system of three boundary 
value problems, defined and coupled over 3 different domains. Moreover, as we 
will see below, the system does not have uniqueness of solution except if the 
parameters involved in the system satisfy a compatibility condition (which is 
satisfied when the parameters come from the full model) and we add an extra 
boundary condition (which is equivalent to set a reference potencial value and, 
therefore, consistent with the physics of the problem). The existence result 
showed in Section 4.3 is not standard, since it requires the use of Functional 
Analysis techniques. 

4.1 Preliminaries 

Let us suppose (ce, Cs, (/ie, </i’s, T) is a solution of the full model. Given t € (0, And) 
let us denote u{x) = (j)eix,t) and v{x) = (j)s{x,t). Then, if 

E{x,w) =exp(^^^{w- f{x))^ -exp^-^^(w-/(a:))^ , 
u and V must satisfy the following system of equations equivalent to (|4|)~([5|): 


{ k-b-{x)E{x, v{x) — u{x)), if a; G (0, Li), 
k^b+{x)E{x, v{x) — u(x)), if X € (Li + S, L), 
0, if X £ (Li, Li + (5), 

ax ax 


(14) 


^(x) =-q-b_{x)E{x,v{x) - u{x)), if x G (0,Li), 
dx^ 

— -^{x)=—q+b+{x)E{x,v{x)—u{x)), if x G (Li + (5, L), (15) 

dx^ 

di’ , du ,, , du , ^ , du , ^ 

— ■^(0) — d-, — “d+, — (Li) — — {Li + (5) — 0, 

ox ax dx dx 


r{x) = £e{xYK{Ce{x,t),T(t)), c(x) = Ce(x,t), V = {I - 2t\) 


o.mt) 


+ > p ’ 


3es(a;) 


b±{,x) = ' {ceix,t)YYcs,±,ma.^ - Cs{x; Rs,t)YYcsix; Rs,t)Y% 

Rs(x) 

q± = 7 -^, d± = /(x) = U{x,Cs{x-,Rs,t),T{t)). (16) 




£s.±CT± A 
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Remark 6 Applying the divergence theorem it is easy to show that any solution 
of [W-CB) must satisfy 

pLi pi 

k- / b-{x)E{x,v{x) — u{x))dx + k+ / b+{x)E{x,v{x) — u{x))dx = 0, 

•^0 Lx-\-S 

pL\ 

q- / b-{x)E{x,v{x) — u{x))dx = d-, 

Jo 

fL 

q+ / b+{x)E{x,v{x) — u{x))dx =—d+. 

Jli+s 

Therefore, a necessary compatibility condition for the existence of solution of 
system is which is true if HM holds, because = 

A ■ 

Remark 7 If we use boundary conditions 111 01) or m the corresponding com¬ 
patibility condition would be or equivalently I ft) = —I ft), which 

is only satisfied if I ft) =0. 

Since function r(x) may not be smooth (at least in a general situation), 
system or, equivalently, ®“®) may not have a classical solution. 

Therefore, we look for what is commonly known as a weak solution. In order 
to find the requirements that this solution shoud satisfy, we use the following 
functional analysis framework. 

Let us consider the Hilbert space 

H = H\0,L) X (^H\0,Li)nH\Li-^6,L)y 


with the norm given by 


(u,v) 


1?^ 


ll^llffqo.L) + ll^llffi( 0 ,Li) + lkllffi(Li+5,L) 




dx 


pLi 

. 2 f dv , .\^ 

/■' 

L 

vO) +(3jW) 

dx + / 

J Z/i -|-5 



dx. 


Here 

€ L^{a, &)| , 

where LP{a, b) (with p > 1) is the set of all measurable functions from (a, b) 
to IR whose absolute value raised to the p-th power has finite integral and ^ 
denotes the derivative of (f) in the sense of distributions. 


H^{a, b) = -Icj) € LJ{a, b) : ^ 
' dx 
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Multiplying (fni) (duality product < >(hi(o,l))'xh^(o,l)) by a function 

(fi e ff^(0, L) and similarly in (ITSl) by a function '0 G Li) fl H^{Li + 6, L) 

we get 








'Li+S 




-k- 

+q- 


b- {x)E(x, v{x)—u{x))ip{x)dx—k^ 


b-(x)E{x, v{x)—u{x))'il^{x)dx+q+ 


b+{x)E{x, v(x) 


JLi+5 

f 

'Li+S 


b+{x)E{x, v(x)- 


-u(x))(p(x)dx 

u(x))ip(x)dx 


= V 



r(x) dc 
c(x) dx 


(x)^(x)dx 

dx 


(i_'0(O) + 


This equality can be rewritten as 


a[(u,v) ,(ip,-ip)} = l[(ip,'ip)j, 

where a [(u, v) , ((p, ip)] and I [(</?, 0)] are the left hand side and right hand side, 
respectively, of the previous expression. According to this, we have the following 
definition of a weak solution (that we will call solution in the following). 

Definition 8 A (weak) solution of system is a couple of functions 

(u, v) € E[ such that 

a [(u, v) , (</?, 0)] = I [((p, 0)] , V {(fi, 0) e H. (17) 


4.2 About the uniqueness of solution of systems (1141) — (I15|) 

For simplicity let us study system (Hl-dlSI), which is equivalent to system (11- 

0 - 

Proposition 9 If {ui,vi), (u 2 ,V 2 } are two solutions of then there 

exists a constant s S IR such that 

U 2 {x) — ui(x) = s, V X € (0, L). 


and 

t> 2 (x) — ui(x) = s, V X G (0, Li) n (Li + 6,L). 

Proof: Let us suppose (ui,ui), {u 2 ,V 2 ) are two solutions of (im) - (fT^ . i.e., they 
satisfy dn. Then, using (</5,0) = {u 2 — ui,V 2 — vi) in (flTt . with 

{ k_ 

— Vi{x) ifxG(0, Li), 

(* = 1 , 2 ), 

—Vi{x) if X G (Li + (5, L), 

q+ 
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we have 


a [( m 2 , ^^2) , (W 2 -Ui,V2- Wl)] - a [(ui, l^l) , {U2 - Ml, V2 - Ml)] = 0 , 
which is equivalent to 


r{x) 


d(M 2 - Ml) 


dx 
2 


(x)) dx 


9- 7o 


Jo V dx J q+ Jl^+s V dx J 


J b-(x)(^E(x,V2(x)-U2(x)) 


-E(x,vi(x) - Mi(x)) j ( (M2(a;) - M2 (x)) - (mi(x) - mi(x)) jdx 


+ fc+ / b+(x)(E(x,V2(x)—U2(x)) 

Jli+S ^ 

— E(x,Vl(x) — Ml(x 


c) — M2 (x)) — (mi(x) — Mi(x)) jdx = 0 . 


Hence, from the monotonicity properties of the exponential function, 


'( 3 :) 


9- Jo V dx J q+ V dx J 

and 

(v 2 ix) — M 2 (x)) — (mi(x) — Mi(x)^ =0 V X S (0, Li) n (Li + (5, L), 
which implies that there exists a constant s G IR such that 
M2(x) — mi(x) = s, V X G (0, L) 


and 

M2(x) — mi(x) = s, V X G ( 0 , Li) n (Li + 5,L). 

Corollary 10 If {ui,vi), (m2,M2) are two solutions of ([1 J| j-| l75|) . then mi(x) — 
Mi(x) = M2(a;) — M2 (x) /or all x G ( 0 , Li) H (Li + d, L) and, therefore, 

E{x, Mi(x) — Mi(x)) = E{x, V 2 {x) — M2(x)). 

Corollary 11 // we add the boundary condition 

m ( 0)=0 ( 18 ) 

(or v(x) = s G IR or m(x) = s G IR, with x and s arbitrarily chosen), then if 
system [Tfl l, (03), (Oil) has a solution, it is unique. 
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Remark 12 According to the mathematical results showed in Proposition 
Corollary [T0| and Corollary [771 a way of getting the property of uniqueness of 
solution of system (gP-0j (see non-uniqueness results in Remark\^ is setting 
a reference potential value for <f>s or (pe for each t S (0,iend)- For instance 

0s(O,i) = 0 (or (psix^t) = s(t) S ]R,x and s{t) arbitrarily chosen). (19) 

4.3 About the existence of solution of system (I14|) . (1151) . 

OS)) 

Let us study the existence of solution for the following linearized approximation 
(see Remark [1]) of system (fill) . (IT5|) . (ITSl) : 


dx 


r(x)^(x) ) + ( r(a:)£ ln(c(a:)) 


dx 


dx 


A_5_(x) (x(x) — u{x) — /(x)) , ii x € (0, Li), 

= { A+6+(a:) (v(x) - u{x) - f{x)) , if x € (Li + 5,L), 
0, if X G (-^1, Li + 5), 


>) = ^(L)=0, 


( 20 ) 


-^(a:) =-6'-6_(x) (■u(x) - u{x) - /(x)), x e (0, Li), 


-^(a:) =-9+b+{x) (v(x) - u(x) - f(x)), x e (Li + 5, L), 






where 


( 21 ) 


;(0) = 0 (or v(x) = s € IR, with x and s arbitrarily chosen), (22) 

, (tta + ac)F (oa + ac)F 

A+ = fc+-, 0+ = q+-. 

Remark 13 The results in Proposition are also valid for the solutions of 
Therefore, if has a solution, then it is unique. 


Multiplypling in (j20p (duality product < •, • >{h^{o,l)Yxh^{o.l)) by a func¬ 
tion (p G H^{0,L) and similarly in the first and second equations of (ITTl) by 

functions -p—'p and respectively, with p G R^(0, Li) fl H^{Li -\- 5,L), the 

U— C7_(_ 

following equality is obtained: 

du d(^ A_ dv dp bw d-)/', , , 
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+A_ 


b-{x) (w(x) — u{x)) {ip{x) — ^(x)) dx 


10 

pL 


+A+ / b+{x) {v{x) — u{x)) {^piyX) — <p(a;)) dx 


/Li+5 


= A_ 


b-(x)f{x) ((p(x) — 'ip(x)) dx + A_| 


/Li+5 


b+{x)f{x) {(f{x) -'ip{x))dx 


, r(x) dc , .do?, ,, X-d- , Aj_fi+ 

+ty I -^ — {x)^{x)dx -2—V'(O) +-5— 


/q c{x) dx ^ ^ dx 
This equality can be rewritten as 


0_ 


a[{u,v) , ((^,V')] = I [(‘/S,-!/')], 

where d [(it, v) , {ip, "0)] and I [{ip, ip)] are the left hand side and right hand side, 
respectively, of the previous expression. According to this, we have the following 
definition of a weak solution (that we will call solution in the following). 


Definition 14 A (weak) solution of system is a couple of functions 

{u, v) € H such that 

a [(u, v) , {ip. Ip)] = I [{ip. Ip)] , V {ip. Ip) e H. (23) 


It is easy to show that, in this case, d : HxH IR is bilinear and continuous 
and I : iJ — > IR is linear and continuous. In order to be able to apply Lax- 
Milgram Theorem we need to show that d : H x H ^ ]Ris coercive, but this is 
not true (otherwise we would get uniqueness of solution of which, as 

we pointed out in Remark 1131 is not true). Let us try to overcome this problem 
by changing iL by a suitable closed subspace. 

Let H = {ip € H^{0,L) : ip{x)dx = 0} x {H^{0, Li) n H^{Li + 6,1)). It 

is known that {ip G H^{0,L) : jJ" ip{x)dx = 0} is a closed subspace of H^{d, 1) 
and there exists a constant a > 0 such that, for any function ip in that set, 

J ^( 2 ^) dx > a J ip{x)‘^dx. 

Proposition 15 There exists a unique couple of functions {u,v) € H, satisfy¬ 
ing 

d [(m, v) , {ip. Ip)] = I [{ip, ip)] , V {ip. Ip) € H. (24) 

Proof: It is easy to show that a : iL x iJ —>■ IR is bilinear and continuous 
and 7 : iL —>■ IR is linear and continuous. Hence, in order to be able to apply 
Lax-Milgram Theorem we need to show that a : ff x iL ^ IR is coercive. Now, 
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nLi i>L 

+A_ / 6_(a:) ('0(x) — (^(x))^ dx + A+ / 6+(x) (■!/'(x) — (/^(x))^ dx 

«/ 0 Z;l+(5 

+ ^ J dx + X-b_J {tf}{x)'^ + ip{x f - 2ip{x)ip{x)) dx 


+A+ 6 _,_ f {tpix)"^ + (p{x)'^ — 2(p{x)ip{x)) dx, 

J L\-\-5 

with r = inf {r(x)} > 0 and = inf { 6 ±(x)} > 0. Using Young inequal- 

xe(0,l) a:G(0,l) 

ity, for any /3_ > 0 and /3+ > 0 we get 


{ip{xf + ip{xf 


2ip{x)ip{x)) > (ip{x) 


I3±il}{xf + ip{xf 



Then, if we take /3_ and /3+ defined by 


Ly*- \ 1 \ 7 


we have that /3_ G ( 0 , 1 ), /3+ G ( 0 , 1 ) and 

a [(¥>, V”) , (V?, V')] > I ^ 

^TrCity) 

-l-A_6_(l —/3_) f ?/;(x)^dx-I-A+6_|_(l —/3+) f ipi^xYda 
Jo Jl,+s 


A_ A 4 


> 10 > Y’ ~ 

which proves the coercivity property and finishes the proof. 


Theorem 16 If (u, v) £ H is the solution of (2^, then {u + ^,v + ^) £ H is a 
solution of \20\) -i21 \} . for any 7 G H. 


Proof: We have to prove that {u, v) = (u-l- 7 , x-l- 7 ) satisfies (l23)) for any 7 G H. 
Now, for any (</?, tp) £ H, 





(p(x)dx. Ip 


1 

L 


ip{x)dx 


£ H. 
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>^-d- _ X+d+ _ I{t) 
A ’ 


Then, using (IMl) and the fact that 
al(u + j,v + j),(<p,^)] 

{u,v),(ip-^J ip{x)dx,'ilj - ^ J ip{x)dx 


= a 


= V5(a;)da:, V' - t / 




X-d- X+d+\ 1 / X , ji ,N 

- J ip{x)dx = I {ip,'ip) 


e_ e+ 

which concludes the proof. 

Corollary 17 There exists a unique solution of r^) - i2S\) . 


5 Conclusions 

A fully mathematical model for the simulation of a Lithium-ion battery is pre¬ 
sented, based on a macro-homogeneous approach developed by Newman (see 
[8]) It includes a system of boundary value problems for the conservation of 
Lithium and conservation of charge in the solid and electrolyte phases, together 
with an initial value problem accounting for the conservation of energy. These 
model can be very helpful for the design and optimization of new batteries and 
also for the real time control of its performance. 

We point out that, over the last years, several authors have used similar 
models with an incorrect boundary condition. In this article it is showed, from 
a physical and a mathematical point of view, why such a boundary condition is 
incorrect. To show that, we prove that the system of equations including that 
boundary conditions does not have any solution. Then, we deduce the correct 
boundary condition and we show some results regarding the uniquesness and 
existence of solution of the corresponding problem. 
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